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Abstract 

A semi-infinite material under a solidification process with the Solomon-Wilson- Alex- 
iades’s mushy zone model with a heat flux condition at the fixed boundary is considered. 

The associated free boundary problem is overspecified through a convective boundary 
condition with the aim of the simultaneous determination of the temperature, the two 
free boundaries of the mushy zone and one thermal coefficient among the latent heat 
by unit mass, the thermal conductivity, the mass density, the specific heat and the two 
coefficients that characterize the mushy zone. Bulk temperature and coefficients which 
characterize the heat flux and the heat transfer at the boundary are assumed to be de¬ 
termined experimentally. Explicit formulae for the unknowns are given for the resulting 
six phase-change problems, beside necessary and sufficient conditions on data in order to 
obtain them. In addition, relationship between the phase-change process solved in this 
paper with an analogous process overspecified by a temperature boundary condition is 
presented, and this second problem is solved by considering a large heat transfer coeffi¬ 
cient at the boundary in the problem with the convective boundary condition. Formulae 
for the unknown thermal coefficients corresponding to both problems are summarized in 
two tables. 

Keywords: Phase Change, Convective Condition, Lame-Clapeyron-Stefan Problem, Mushy 
Zone, Solomon-Wilson-Alexiades Model, Unknown Thermal Coefficients. 

2010 AMS subjet classification: 35R35 - 35C06 - 80A22 


1 Introduction 

Heat transfer problems with a phase-change such as melting and freezing have been studied in 
the last century due to their wide scientific and technological applications. Some books in the 
subject are [HIHj6l H0lfl2lH6l 1261135] . 

In this paper we consider a phase-change process for a semi-infinite material, which is char¬ 
acterized by x > 0, that is initially assumed to be liquid at its melting temperature (which 
without loss of generality we assume equal to 0°C). We consider this material under a solidi¬ 
fication process with the presence of a zone where solid and liquid coexist, known as ’’mushy 
zone”, with a heat flux boundary condition imposed at the fixed face x — 0. We follow [301133] 
in considering three different regions in this type of solidification process: 
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1. liquid region at temperature T(x,t ) = 0: A = {(a;, t) G M 2 / x > r(t), t > 0}, 

2. solid region at temperature T(x,t ) <0: D s = {(x,t) G M 2 /0 < x < s(t), t > 0}, 

3. mushy region at temperature T(x,t ) = 0: D p = {(x,t) G M 2 / s(t) < x < r(t), t > 0}, 

being x = s(t) and x = r(t) the functions that characterize the free boundaries of the mushy 
zone. We also follow [30] in making the following assumptions on the structure of the mushy 
zone, which is considered as isothermal: 

1 . the material contains a fixed portion of the total latent heat per unit mass (see condition 
© below), 

2 . its width is inversely proportional to the gradient of temperature (see condition (ED below). 

Parameters involved in this problem are: 

l > 0 : latent heat by unit mass, 
k > 0 : thermal conductivity, 
p > 0 : mass density, 
c > 0 : specific heat, 

0 < e < 1 : one of the two coefficients which characterize the mushy zone, 

7 > 0 : one of the two coefficients which characterize the mushy zone, 
q 0 > 0 : coefficient that characterizes the heat flux at x — 0 , 
ho > 0 : coefficient that characterizes the heat transfer at x — 0 , 

—D 0 o < 0: bulk temperature at x — 0. 

We suppose that five of the six thermal coefficients /, k, p, c, e and 7 of the solid phase are 
known and that, by means of a change of phase experiment (solidification of the material at its 
melting temperature) we are able to measure the quantities qo, ho and — Ax>- 

Encouraged by the recent works [3T1134] and with the aim of the simultaneous determination 
of the temperature T = T(x, t), the two free boundaries x = r(t ) and x = s(t), and one unknown 
thermal coefficient among l, k, p, c, e and 7 , we impose an overspecified boundary condition [4] 
which consists of the specification of a convective condition at the fixed face x = 0 (see condition 
© below) of the material undergoing the phase-change process. This lead us to the following 
free boundary problem: 


pcT t (x,t ) - kT xx (x,t) = 0 

0 < x < s(t), t > 0 

( 1 ) 

T{s(t),t ) = 0 

t > 0 

( 2 ) 

kT x (s(t),t) = pl[es(t) + (1 - e)r(f)] 

t > 0 

(3) 

T x (s(t),t)(r(t ) - s{t)) = 7 

t > 0 

(4) 

r( 0 ) = s( 0 ) = 0 


(5) 

kT,(0,t) = ^= 

t > 0 

( 6 ) 

M’ x ( 0 ,i) = ^(T( 0 ,«) + C oo ) 

t > 0 

(7) 


This problem was first study in pJ3] with a temperature boundary condition at x = 0 
instead of the convective condition (J7D considered in this paper. Moreover, the determination 
of one unknown thermal coefficient for the one-phase Lame-Clapeyron-Stefan problem with an 
overspecified heat flux condition at the fixed face x = 0 without a mushy zone was done in [32] . 
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Other papers related to determination of thermal coefficients are 0|3llMl[IIllIM5llI7H25l[23- 

mmm. 

The goal of this paper is to obtain the explicit solution to phase-change process (HD-CD) 
with one unknown thermal coefficient, and the necessary and sufficient conditions on data 
in order to obtain an explicit formula for the unknown thermal coefficient. In addition, we 
are interested in analysing the relationship between problem ®-o and the phase-change 
process given by ([I])-© beside the Dirichlet boundary condition overspecified at x = 0 given 
by Q35D (see below). In particular, we are interested in solving the problem with Dirichlet 
boundary condition through problem with convective boiindary condition when large values of 
the coefficient ho that characterizes the heat transfer at x = 0 are considered. 

The organization of the paper is as follows. In Section [2] we prove a preliminary result where 
necessary and sufficient conditions on data for the phase-change process (JT j) -(JT| ) are given in 
order to obtain the temperature T = T(x, t) and the two free boundaries x = r(t) and x = s(t). 
Based on this preliminary result, in Section |3] we present and solve six different cases for the 
phase-change process ([I])-© according the choice of the unknown thermal coefficient among 
l, k, p, c, e and 7 . In Section [4] we discuss the relationship between the phase-change process 
m-m with the Dirichlet boundary condition (135]) and the same process with the convective 
boundary condition ©. We show that the temperature T D = T D (x,t), the free boundaries 
x = and x = S£>(t), and the explicit formula for the unknown thermal coefficient l , k, p, 

c, e or 7 for the phase-change process (fTj) -(j 6 j) with the Dirichlet condition (j35j) can be obtained 
through the phase-change process with convective condition given by (]Tj)-(JTJ) when ho tends to 
+ 00 . Explicit formulae for the unknown thermal coefficient for problems ©-© and o-o 
and (155]) . beside restrictions on data that guarantees their validity, are summarizes in Table 1 
and 2 respectively. 


2 Explicit solution to the phase-change process 

The following lemma represents the base on which the work in this Section will be structured. 
Lemma 2.1. The solution to problem m-m is given by: 

q 0 \Jtux 


T(x, t) = 

Ah 

s(t) = 2^\fcd 


erf ( — erf(£) 


\2 s/at J 


r(t ) = 2 




7 k 


2 q 0 y/a 


exp (£ 2 


y/at 


0 < x < s(t), t > 0 

( 8 ) 

t > 0 

(9) 

t > 0 

( 10 ) 


if and only if the parameters involved in problem ©)-([7|) satisfy the following two equations: 


£ + 


7^(1 — e) 


2q oy /a 
kD < 

erf(f) = 


exp (£ 2 


exp (£ 2 ) = 


% 


ply/a 


% 


qoy/na 


where a = ^ represents the thermal diffusivity. 


hr\D 


O-^oo 


( 11 ) 

( 12 ) 


Proof. The kind of phase-change processes considered in this article have the following general 
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solution [801I8T1I88] : 


n*,t)- a+b*^) 

0 < x < s(t),t > 0 

(13) 

s(t ) = 2 £\/af 

t > 0 

(14) 

r(t) = 2 p\fcrt 

t > 0 

(15) 


where coefficients A, B, £ and p depend on the particular specifications of the phase-change 
process. 

I11 orden to have the solution to problem ©-(IT]), we impose conditions (J2])-(J4J) , ([6]) and (Hj) 
on ffT3]) -(IT5 ]) and obtain that coefficients A, B and p must be given by: 


A 



erf(O) 


B 


qoy/na 

k 


and 


yfcexp (£ 2 ) 
2 qo^fa 


(16) 


which corresponds to solution flHD-OQD, and the parameters involved in the problem must satisfy 
equations (Hill and (TT2jh ■ 

As a consequence of Lemma 12.11 we know that we can solve the phase-change process 
(II])-© with one unknown thermal coefficient through the determination of the parameter £ 
that characterizes one of the two free boundaries of the mushy zone and the unknown thermal 
coefficient among l, k, p , c, e and 7. In addition, we also know from Lemma 12.11 that we can 
do that by solving the system of equations (TTOl - ffTSl) . 


3 Explicit formula for the unknown thermal coefficient 

In this Section we present and solve six different cases for the phase-change process (TTT)- (T7T) 
according the choice of the unknown thermal coefficient among l, k, p, c, e and 7. 

With the aim of organizing our work, we classify each case by making reference to the co¬ 
efficients which is necessary to know in order to solve it (see Lemma EZQ): 

Case 1: Determination of l and £, Case 2: Determination of 7 and £, 

Case 3: Determination of e and £, Case 4: Determination of k and £, 

Case 5: Determination of p and £, Case 6: Determination of c and £. 


In addition, with the goal of make our presentation more readable, in the following statements 
and proofs we introduce several functions. We name this functions with a subscript according 
the case where they arise. 

Theorem 3.1 (Case 1: determination of / and £). If in problem (QJ)-([?p we consider the thermal 
parameter l as an unknown, then its solution is given by /f^l)-/f7dl) with l and £ given by: 



(17) 

(18) 
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if and only if the parameters qo, h 0 , D aQ , k, p and c satisfy the following two inequalities: 

% 


1 - 


hnD 


> 0 


Doc 

do 


O-^oo 

kpc 


71 


1 - 


% 


hnD 


O-^oo 


< 1. 


(Rl) 


(R2) 


Proof. Due to properties of the error function, it follows that a necessary and sufficient condition 
for the existence and uniqueness of a positive solution to equation (fT 2 ]l is: 


0 < 


D r 


do 


kpc 

71 


1 - 


% 


hnD 


O-^oo 


< 1, 


which is equivalent to inequalities (IRip and 


In that case, the positive solution to equation 


(fT2j) is given by (1 1 81) . Finally, it follows from equation (flip that l is the positive thermal 
coefficient given by (IT7|) . ■ 

Theorem 3.2 (Case 2: determination of 7 and £). If in problem J7])-([?]) we consider the thermal 
parameter 7 as an unknown, then its solution is given by m-m with 7 given by: 


7 = 


2 go 


go 


/ 2 (0 exp (— 2 £ 2 ) 


(l-e)VWc\l\ Pk JZVS "^ V ~ S ' (19) 

and f given by H8\) . if and only if the parameters qo, ho, -Doo, k, p, c and l satisfy inequalities 
\R1\) . and: 

'Dr 


h erf 


-1 


kpc 


go 

hoDoo 


< 


go 

l 



c 

pk’ 


(R3) 


9o * IT 

where the real function / 2 is defined by: 

/ 2 (x) = a: exp (x 2 ), x > 0 . ( 20 ) 

Proof. As we see in the proof of Theorem 13.11 a necessary and sufficient condition that guar¬ 
antees the existence and uniqueness of solution to equation (TT2|) is that inequalities (1Rl[i and 
(IR.2D hold, and in that case, the coefficient f is given by 

On the other hand, it follows from equation (ITT]) that 7 is given by (fTOlh This coefficient is 
positive if and only if: 

< T^k’ (21) 

where / 2 is the real function defined in (1201) . Taking into account the expression of f given in 
(in, we have that inequality (ITT]) is equivalent to inequality 


Theorem 3.3 (Case 3: determination of e and £). If in problem (Q]j-([^) we consider the thermal 
parameter e as an unknown, then its solution is given by m-m with e given by: 


e = 1 


2 go 


go 


7 \Jkpc V l 



/ 2 (f) exp(- 2£ 3 


( 22 ) 


and f given by HM . if and only if the parameters qo, ho, D 0Q , k, p, c and 7 satisfy inequalities 
ITri\) . nm) . mm and: 


7 Vkpc 


2 go 


exp 


/2 erf 


-1 


erf 

D 


-1 


Dr 


go 


kpc 

71 


go 


hoD, 


O-^oo 


go 


kpc 

71 


1 - 



(R4) 
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where fo is the real function defined in K20 j) . 

Proof. Conditions (IR.lj) and (1R.2[) . and the expression of £ given in (1T81) arise in the same way 
that in the precedent proofs. On the other hand, it follows from equation (1111) that e is given 
by (l22j) . being / 2 the real function defined in (1201) . This coefficient is positive if and only if: 


2go 

"fVkpc 


ffo 

\ i y pk 



exp (—2£ 2 ) < 1. 


(23) 


Taking into account the expression of £ given in (ITS]) , we have that inequality (1251) is equivalent 
to inequality (IR dj) . Finally, we have that e given in (1221) is less than 1 if and only if /2(£) < , 

which, as we see in the proof of Theorem 13.21 is equivalent to condition (I R3 1) . ■ 


Theorem 3.4 (Case 4: determination of k and £). If in problem we consider the thermal 

parameter k as an unknown, then its solution is given by m-m with k given by: 

r -i 2 


k 


7T 

pc 


9oerf({) 



(24) 


and £ the unique solution of the equation: 



where / 4 is the real function defined by: 


% \ 

hoDoo J 


x > 0, 


(25) 


U{x) = 


X + 


70 r(l 


2 D c 


<20 


HqDc> 


erf (a:) exp (x 2 


erf(x) exp (x 2 ), x > 0, 


(26) 


if and only if the parameters q 0 , h 0 and satisfy inequality m- 
Proof. The system of equations (fTTj) - (fT2j) is equivalent to: 

(27) 


(28) 

where / 4 is the real function defined in (126|) . A necessary condition for existence of solution to 
this system is that inequality (IR1I) holds. Then, if we assume that (IR 1 jl holds, we inmediately 
obtain that k is given by (12T|) . To complete the proof only remains to demonstrate that equation 
(125|) admits a unique positive solution. This follows from the fact that / 4 is an increasing 
function such that /(0 + ) = 0 and /(+oo) = +oo. ■ 


\fk = 

m) - 



g 0 erf(£) 


P c D, 

C-D oo 

ly/n 


1 - 


<?0 


hoD 0 

Qo 


hnD , 


0 


Theorem 3.5 (Case 5: determination of p and £). If in problem (QJ)-([7p we consider the thermal 
parameter p as an unknown, then its solution is given by Iffll-KTU) with p given by: 


P 


7T 

kc 


Qo er f(£) 


D oo ( 1 


90 


h()Do 


(29) 


and £ the unique solution of the equation (H3J), if and only if the parameters qo, h 0 and D^ 
satisfy inequality m- 
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Proof. It is similar to the proof of Theorem 13.41 


Theorem 3.6 (Case 6 : determination of c and £). If in problem |IJ)-([?p we consider the thermal 
parameter c as an unknown, then its solution is given by m-m with c given by: 

r i 2 


c 


7T 

pk 


®erf({) 



(30) 


and f the unique solution of the equation: 


h(x) = 


QoV^ 


plkD, 


1 _ go 

00 ' h 0 D o 


x > 0, 


(31) 


where f e is the real function defined by: 


h(x) = 


x 


70r(l -e) 


erf(x) 


2D^ 1 


go 


hoD 0 


exp (x 2 ) 


exp (x 2 ), x > 0 , 


(32) 


if and only if the parameters qo, h 0 and D 0a satisfy inequalities m and: 


i-^-< 1 




- 7 ( 1 - 6 ) 


^ 0 -^CXD -^00 _plk 

Proof. The system of equations ffTT]) - ffT 2 ]) is equivalent to: 


(R5) 


V~C=\J 

r*~ Qo erf (f) 

(33) 

pk d (l 90 ^ 

U 00 V hoDoc) 

h{x) = 

qIV^ 

(34) 


plkDce (l hotf) 


where / 6 is the real function defined in (I32p . A necessary condition for existence of solution to 
this system is that inequality (1R 1 1) holds. Then, if we assume that OR 111 holds, we immediately 
obtain that c is given by (1301) . To complete the proof only remains to demonstrate that equation 
(13T)) admits a unique positive solution. Since fe is an increasing function such that /( 0 + ) = 
f + \ and /(+ 00 ) = + 00 , it follows that a necessary and sufficient condition for 

2D °° V 1_ hQDaa ) 

existence (and uniqueness) of solution to equation (l3lj) is that: 

7 T 7(1 -e) 0 F < q 2 ^ 

20)00 (l - plkDoo (l - 

which is equivalent to inequality (IR5I) . ■ 

Table 1 summarizes the results of this section, corresponding to 6 cases. 
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Case 


Thermal coefficient 


Coefficient £ that characterizes 
the free boundary x — s(t) 


Restrictions on data 


/ rr go exp (-g 2 ) 

- Y pk ^ + 7 (l- 0 ^ exp( g 2 )j 


£ = err 1 (i- 

^ ' 90 V 77 V 


90 


ho Do 


m, w 


_ 2 g 0 / go 


^ ( 1 — e)-\/kpc l l y pk 


— £ exp (£ 2 ) ) exp (— 2 £ 


£ = erf ” 1 1 - 

s ' 90 V n V 


90 


ho Do 


(ED, (E2D, (11 


e = 1 - f a/ 5 - ^ ex p (a) ex p (- 2 a 


£ = erf ” 1 1 

^ ' 90 V ^ \ 


go 


hoD a 


m, W, (ESD, dH 


Unique positive solution of: 


k = — 

pc 


90 erf (Q 


Doo 1 - 


HqDq 


h{x) = 


1 - 


lyj 7T \ h 0 D 


, with / 4 defined by 


m 


Unique positive solution of: 


P = 


kc 


go erf(£) 


h( x ) = ( 1 - TUzb)’ with / 4 defined by (® 


m 


c =^ 


Qo er f (U 


£>oo 1- 


90 


ho D 0 


Unique positive solution of: 

M x ) = ,, n if* 90 \ > with h defined by 

plkU oo^i hQDoo ) 


m, (ESD 


Table 1: Formulae FOR PROBLEM ©-©■ Explicit formulae for the unknown thermal coefficient l, 7 , e, k, p or c and coefficient £ 
(or the equation that it must satisfy) and the corresponding restrictions on data that guarantees their validity. 





































4 The phase-change process with large heat transfer co¬ 
efficient 


A similar phase-change process to o-o with one unknown thermal coefficient have been 
studied in [33]. In that paper, the author consider a fusion process with a mushy zone given by 
©-© overspecified with a temperature boundary condition and obtain the explicit solution 
to some cases. Encouraged by [33], let us consider the solidification process CO)-© with one 
unknown thermal coefficient overspecified with the Dirichlet boundary condition: 


T( 0 ,f) = -.Doo, t > 0 . 


(35) 


We can see this condition as the limit case of the convective boundary condition © when 
the heat transfer coefficient ho tends to +oo. From a physical point of view, overspecify the 
phase-change process ©-a by imposing the convective boundary condition © seems to be 
more appropriate than imposing the Dirichlet condition (135 D . This Section is devoted to show 
that the temperature T D = T D (x,t), the free boundaries x = r D {t ) and x = sn(t), and the 
explicit formula for the unknown thermal coefficient l , k , p, c, e or 7 for the phase-change 
process with Dirichlet boundary condition given by ©-© and (1351) can be obtained through 
the phase-change process with convective boundary condition given by ©-© when ho tends 
to + 00 . 

We begin with a result related to the solution to the phase-change process ©-© and (I33|) . 
This result may be shown in much the same manner as Lemma 12.11 thus we do not give its 
proof here. 


Lemma 4.1. The solution to problem ©)-©) and TTR) is given 

<Wtto 


T D (x,t ) = 


k 


erf ( 2 ^) “ erfK) 


0 


su(t) = 2 fy/at 

r D (t) = 2 


c+A^ ex p« 2 ) 

2 qoy/a 


y/od 


by nto-HTh). that is: 

< x < sr>(t), t > 0 
t > 0 
t > 0 , 


where a = ^ represents the thermal diffusivity, if and only if the parameters involved in problem 
©)-©) and ( 0 © satisfy the following two equations: 


? + 7(1 2r , £)7i: erf (f) exp K 2 ) 

erf( { ) = E?L.[E*. 

qo V 7T 


erf (0 exp (e 2 ) 


C-Dqo 


W vr 


(36) 

(37) 


From Lemma 14.11 we know that, in order to have the temperature Tp = Tjj(x,t), the free 
boundaries x = rn[t) and x = s_o(t), and the unknown thermal coefficient l, k, p, c, e or 
7 for problem ©-© and ([33]) . it is enough to find the unknown thermal coefficient and the 
parameter that characterizes the free boundary sn(t). Proceeding analogously to the work done 
in [33] or in Section [3] we can obtain the thermal coefficient l, k, p, c, e or 7 and the parameter 
£ for problem ©-© and f[35l) . Formulae for those quantities, besides restrictions on data that 
guarantee their validity, are summarized in Table 2 (restrictions on data and definitions on 
functions mentioned in Table 2 are listed below the table). 
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Case 

Thermal coefficient 

Coefficient £ that characterizes 
the free boundary x = sn(t) 

Restrictions on data 

1 

? rr <?o exp (-£ 2 ) 

£ = erf -1 ( 
s \ 90 V ^ y 

m 

- y pk [g+^-^exp^)] 

2 

— 2 9o {90 rr 

-eexp (e 2 )) exp (- 2 a 

f = erf- 1 ( ./Hp£\ 

^ \ Q0 \ TV 1 

(IR7D 

' (1 —e)y/kpc \ l V pk 

3 

' = 1 - vis (ly^ " 5 exp « 2 )) “P (-25 2 ) 

£ = erf -1 ( 

^ \ qo \ tv 1 

(lR7l). (lR8l) 

4 

k = — 

pc 

<?oerf(£) 

Doo 

2 

Unique positive solution of: 

F 4 (x) = with F 4 dehned by fl40l) 


5 

p=fc 

90 erf(£) 

Dqg 

2 

Unique positive solution of: 

F a {x) = with F 4 dehned by fl40l) 


6 

C = ~pk 

90 ef f(?) 

Doo 

2 

Unique positive solution of: 

f 6 {x) = Wlth Fe defined b y © 

f[R9l) 


Table 2: FORMULAE FOR PROBLEM (fTj) -(l6l) AND (1351) . Explicit formulae for the unknown thermal coefficient l, 7, e, k, p or c 
and coefficient £ (or the equation that it must satisfy) and the corresponding restrictions on data that guarantees their validity. 
































List of restrictions on data for problem (TTj) -(]6| ) and fl35l) mentioned in Table 2: 


D c 


% 
D « 


kpc 


7T 


kpc 


< 1 


< erf (rj) 


qo V 7r 

where rj is the unique positive solution to the equation: 


f i \ Qo c 

h(x)= -q-pk' x>0 ' 


being / 2 the real function defined in (J2UD - 


(R6) 

(R7) 


(38) 


kpc 


> erf (ri), x > 0, 

?0 * IT 

where q is the unique positive solution to the equation: 


(R8) 


h(x) 


^y/kpc 2 q 0 \~C 
— exp (2*) = -^, 

being / 2 the real function dehned in (l20|l . 


x > 0, 


(39) 


lkpD c 
2 qo 


1 + 


7(1 - e) 


< 1 


(R9) 


List of definitions of functions related to problem (fT jl -lfh ]) and (1HHD mentioned in Table 2: 


F 4 (x) = erf (x)f 2 (x) + 


(1 - 6)7 
2D^ 


[erf (x)] 2 exp (2x 2 ), x > 0, 


(40) 


being / 2 the real function dehned in (f2Ul) . 


FcM= (sm + 7 ( 1 2 c!! v? exp (j2) ) exp (j2) - i>0 (41) 

On the other hand, it is not difficult to verify that formulae and restrictions on data given in 
Table 2 correspond to formulae and restrictions on data given in Table 1 for ho tending to + 00 . 
This fact, besides Lemmas 12.11 and 14.11 allow us to conclude that we can solve the phase-change 
process 0-0 with one unknown thermal coefficient overspecified by the Dirichlet condition 
(j55|) through the phase-change process (JT j) - ()TI) . which is overspecified by the more physically 
appropriate convective boundary condition (]7|). when the heat transfer coefficient h 0 tends to 
Too. 


5 Conclusions 

In this paper, we consider a semi-infinite material under a solidification process with a mushy 
zone caused by an initial heat flux boundary condition. We solve the associated free boundary 
problem overspecified with a convective boundary condition and obtain the temperature, the 
two free boundaries of the mushy zone and one thermal coefficient among the latent heat by 
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unit mass, the thermal conductivity, the mass density, the specific heat and the two coefficients 
that characterize the mushy zone, when the bulk temperature and the coefficients that char¬ 
acterize the heat flux and the heat transfer at the boundary are assumed to be known. As a 
consequence, we give formulae for the temperature, the two free boundaries and the unknown 
thermal coefficient, beside necessary and sufficient conditions on data in order to obtain them. 
In addition, we present the relationship between the phase-change process studied in this pa¬ 
per with another similar phase-change process which is overspecified by a Dirichelt boundary 
condition. From this relationship, we solve the problem with the Dirichlet condition by con¬ 
sidering a large heat transfer coefficient in the problem with the convective condition. In this 
way, we solve the phase-change process overspecified with a temperature boundary condition 
through the more physically appropriate phase-change problem overspecified with a convective 
boundary condition. We summarize explicit formulae for the unknown thermal coefficient for 
both problems in Tables 1 and 2. 
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